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Abstract 

We consider an autonomous differential system in M n with a periodic 
orbit and we give a new method for computing the characteristic multipliers 
associated to it. Our method works when the periodic orbit is given by the 
transversal intersection of n — 1 codimension one hypersurfaces and is an 
alternative to the use of the first order variational equations. We apply it 
to study the stability of the periodic orbits in several examples, including a 
periodic solution found by Steklov studying the rigid body dynamics. 
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1 Introduction and statement of the results. 



Consider a differential system in R n , with n > 2, given by: 

I - X(x), (1) 

where X:WC R n —>■ R™ is a C l function in some non-null open set WC1™ and t 
is a real independent variable. We assume that system (JTJ exhibits a periodic orbit 
T := {^(t) | < t < T} C U with period T > 0. As usual we will denote by 4>(t,p) 
the flow solution of (JTJ) such that 4>(0,p) = p, by Xi the 2* h component of the point 
x, that is, x = (xi, X2, ■ ■ ■ , x n ) T and analogously X = (X±, X 2 , . . . , X n ) T , where T 
denotes transposition. 

It is well known that to determine the behavior of the flow near T, a first step 
is to get the characteristic multipliers associated to this orbit. These multipliers 
are usually obtained through the study of the first order variational equations. In 
this paper we propose an alternative way for obtaining them which works when T 
is given by the transversal intersection of n — 1 codimension one hyper surf aces. 

Recall that the behavior near V is given by the Poincare map, which is defined in 
a section E through a point p G T, where E passes through p and is a local smooth 
manifold of dimension n — 1 transversal to V. Given a section E, the Poincare map 
is defined as: 

II : E — ► E 

q h-* (t>(r(q),q), 

where r{q) is the unique real function such that 0(r(g), q) G E and limr(g) = T. 

The n — 1 eigenvalues of DT\.(p) are independent of p and E and are called the 
characteristic multipliers of Y. It is well known that the stable (respectively, unsta- 
ble) manifold associated to V has dimension the number of characteristic multipliers 
with modulus smaller than 1 (respectively, bigger that 1). It is also well known, 
see jBj, that if all the characteristic multipliers have modulus lower than or equal to 
1, then T is Liapunov stable and if all the characteristic multipliers have modulus 
strictly lower than 1, then it is asymptotically stable. 

Our result is motivated from a previous result given in [H] for planar differential 
systems. Recall that in the planar case there is only one characteristic multiplier 
which is given by: 

T 



n'(p) = ex p|^ div(X)(7(*))tft 



where div(X) = 8X1/8x1 + 8X2/8x2 is the divergence of the system, see for instance 
[TUl p. 214]. In the authors give an alternative expression for the same value 
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n'(p). Let us consider an invariant curve f{x\,X2) = for a planar system 
that is a curve defined by a real C 1 (W)-function f{xx,x^), for which there exists a 
function 22) : W C I 2 -> 1 of class C 1 satisfying that: 

D/(x) X(x) = V/(x) ■ X(x) T = fe(x) /(x), 

where V/ = (df/dx\, df/dx^)- Under these assumptions, the function /c(xi, ^2) is 
called the cof actor associated to the invariant curve f(x\,X2) = 0. 

Proposition 1 [5J Consider system (JIJ) wift n = 2 and a periodic orbit T = 
{ 7 (t) I < £ < T} contained in a planar invariant curve f(xi,X2) = 0. // V/ 
does not vanish on T, then 

T 



II'(p) = expj^ fc( 7 (*))d*}. 



The above result allows to compute the characteristic multiplier associated to T 
through the integration over the cofactor associated to an invariant curve contain- 
ing the periodic orbit. In Theorem El we extend this result to periodic orbits in W 1 
given as the transversal intersection of n — 1 codimension one hypersurfaces. 



Theorem 2 Let V = { 7 (i) | < t < T} be a T— periodic solution of (|T|). Consider 
a smooth function f : W C 1" -> M ra_1 , f = / 2 , . . . , / n _i) T , suc/i i/ioi: 



T is contained in [^j {/i(x) = 0}, 



i=l,2,...,n-l 

- t/ie crossings of all the manifolds {/i(x) = 0} /or z = 1,2, ...,rt — 1 are 
transversal over T, 

- there exists a (n — 1) x (n — 1) matrix k(x) 0/ rea/ functions satisfying: 

Df(x)X(x) = k(x)f(x). (2) 

Let v(i) 6e t/ie (n — 1) x (n — 1) fundamental matrix solution of 

^ = k( 7 (t))v(t), v(0) = ld. (3) 
Taen tae characteristic multipliers ofT are the eigenvalues o/v(T). 
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Indeed in the proof of this theorem, we will deduce that the matrices v(T) and 
DU(p) are conjugated. On the other hand, recall that the usual way of comput- 
ing the characteristic multipliers of F lies in finding the fundamental matrix u(T), 
solution of the variational equation: 

^ = DX(7(t))u(t), u(0) = ld, (4) 
which is a n x n matrix having eigenvalues 1 and the n — 1 characteristic multipliers 

of r. 



Note that the assumptions stated in our theorem are straightforward generaliza- 
tions of the notion of invariant curve containing the periodic orbit V for a system 
in dimension n and the integration of the cofactor, which in R n consists on solving 
the linear differential system v'(£) = k(7(i)) v(i). 

Next section is devoted to prove the above theorem. Finally, in Section |3] we 
apply it to study the stability of the periodic orbits of several differential equations. 
Example 1 is devoted to study a 4-dimensional polynomial system, that includes 
two systems studied in exhibiting an explicit periodic orbit. Example 2 shows a 
3-dimensional polynomial system for which we prove that the stability of the given 
periodic orbit is equivalent to the study of the stability of the Mathieu's equation. 
The last example deals with a more involved 6-dimensional system that 

controls the dynamics of a rigid body. In our approach, and by using the existence 
of three independent first integrals, we reduce the study of the stability of the 
Steklov periodic orbit to the study of a second order linear differential equation. 
This Steklov orbit is introduced in [TT]. It is worth to say that in all the above 
examples we have also tried to use the usual approach, namely the variational 
equations, for studying the stability of the given periodic orbit and we have found 
that the result of Theorem |21 makes the computations easier. 



2 Proof of Theorem [2] 

Fix a point p G T. To prove the result we consider the Poincare section given by 
the (n — 1)— dimensional orthogonal hyperplane to X(p), 

£ =< X(p) >^=< ei, . . . e n _! > . 

Take a new system of coordinates centered at p with basis ei, . . . e n _i, X(p). 
These new coordinates write as y = A(pc — p), for some invertible matrix A. By 
using them the differential equation JI| is converted into 

% = Y(y):= AX(A-iy+p), (5) 
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and the manifold containing T is given by g(y) := f(A 1 y + p) = 0, which satisfies 

£g(y)Y(y) = k(y)g(y), (6) 
where k(y) := k(v4 _1 y + p), because by using ©, the following equality holds: 
Di (A _1 y + p) A' 1 A X(A _1 y +p) = k(A _1 y + p) f (A _1 y + p). 

Recall that <f)(t,p) is the solution of ([!} such that when t = passes through 
p. Let ip(t,q) be the solution of (0) such that when £ = passes through q. Then 
V>(*,g) =A[<f){t,A- 1 q+p)-p]. 

In these new coordinates note that E = {y n = 0} and if g = (yi, y 2 , ■ ■ ■ , y n -i, 0) T € 
E then the Poincare map EE : E — > E writes as 

n(g) = (>i(>(g), g), ^2(^(5), ?),•••, (?"(?)> 9), 0) T , 

where r(g) is precisely the time such that ip n {{T~(q), q)) = 0, which is known to be 
a smooth function. Indeed we can identify II with a map II from R n ~ l into itself 
with variables yx, . . . , y n -±, defined as 

U(y ly . . . ,y n -i) = (ITi(?/i, . . .,y ft _i,0), . . . , IIn_i(j/i, . . . ,y„_i,0)) . 

Observe also that the hypotheses on the transversal cutting of T and the hy- 
persurfaces /i(x) = imply that extended map, f = (fi, f 2 , . . . , f n -i, f n ), where 
/ n (x) := Xi{p)x\ + X 2 {p)x 2 + ■ • ■ X n (p)x n is such that det(£)f(p)) 7^ 0. This infor- 
mation translated to the function g(y) implies that the matrix -Dg(O) has rank n — 1 
and that a minor with determinant different from zero is the one corresponding to 
the partial derivatives with respect to yx, y 2 , ■ ■ ■ , y n -x- 

On the other hand, by using © we have that: 

^<|M> = D B W, g) ) ^fi» = Dg (*(*,,)> Y(#,,)) 

= kW>(f,g)) g «■(<,?)), 

for any point g in the domain of definition of (j3J). 
Let v(t; g) be the fundamental matrix solution of 

^ = k(V(t,g))v(t). (7) 

Thus g (■?/>(£, g)) = v(£; g) g(g), because the function g (if)(t, q)) is a solution of the 
same linear differential system and satisfies that g (^(O, g)) = g (g) . 

Consider a point g £ E. In these coordinates, we have g = (gi, q 2 , ■ ■ ■ , q n -i, 0) T . 
We define z = (gi, g 2 , . . . , g n -i) T and h(z) = (z T , 0) T = g like the inclusion of z as 
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a point in the domain of definition of (j5J). Since H(q) = V'C^C'?) 5 ^) 5 by the above 
result we have that g (11(g)) = v (r(g), q) g(z), or equivalently that 



h(U(z)))=v(r(h(z));h(z)) g(h(z)). 

We differentiate the previous identity with respect to z: 

D S (h(U(z))^ Dh(U(z))Df[(z) = D[v(r(h{z));h(z))} g(h(z)) + 

+ v(r(h(z));h(z)) Dg{h(z))Dh(z). 

By evaluating at z = 0, which corresponds exactly to the point p G T, and using 
that g(h(0)) = 0, /i(0) = (0 T , 0) T , r(0) = T, ft(0) = and the expression of h, we 
obtain that: 

D zg (0)DIL(0) = v(T;/i(0)) D,g(0), 

where D z g(0) := Dg(/i(0)) Dh(0) is precisely the invertible squared matrix formed 
by the derivative of g only with respect to qi, ■ ■ ■ ,q n -i- Hence we have that the 
characteristic multipliers associated to T coincide with the eigenvalues of v(T; 0) 
because the matrices DU(0) and v(T; h(0)) are similar. 

Finally note that k(ip(t, h(0))) = k(A -1 ^(t, h(0))+p) = k(</>(t,p)). Thus equa- 
tion (J7J) is the same that 

and v (T; /i(0)) coincides with the matrix v(T) as defined in (JSJ), as we wanted to 
prove. 



3 Some examples and applications. 

Example 1. The following example is an extension of two systems extracted 
from 4j. The goal of the examples given in 0] is to illustrate that the asymptotic 
stability of a periodic orbit of a system (0) is not determined by the eigenvalues of 
the matrix defining the first variational equation. 

We consider the following differential system in R 4 : 

x = —y — x (x 2 + y 2 — l) , 

y = x - y (x 2 + y 2 - 1) , 

z = —w — s z {z 2 + w 2 — l) — s k (x — z) , 

w = z — sw [z 2 + w 2 — l) — s k (y — w) , 

with s, k real parameters. This system coincides with the first example given in 
0] when s = 1 and with the second example when s = —1. We note that 
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this system always exhibits the periodic orbit T = {"f(t) : < t < 2n} , where 
7(t) = (cos(t), sin(t), cos(t), sin(t)). We are able to compute all the characteristic 
multipliers associated to T for any real value of the parameters s and k. We consider 
the hypersurfaces fi(x, y, z, w) — 0, i — 1, 2, 3, given by: fi(x,y,z,w) = x 2 + y 2 — l, 
f 2 (x,y,z,w) = x - z and f 3 (x,y,z,w) = y-w. We denote by f = (/i,/ 2 ,/ 3 )- 
It is easy to see that T is contained in the intersection of these three hypersurfaces 
and that the crossings of these hypersurfaces are normal over the periodic orbit T, 
as the computation of the following determinant shows: 



4^0. 



V/i(7(*)) 




2 cos(t) 


2 sin (t) 








V/ 2 (7(*)) 




1 





-1 





V/ 3 (7(*)) 







1 





-1 


X( 7 (0) 




— sin(t) 


cos(t) 


— sin(i) 


cos(t) 



Straightforward computations show that Df (x) X(x) 
ing matrix of cofactors: 



k(x) f (x), with the follow- 



/ -2 (x 2 + y 2 ) 



k(x,y,z, w) :-- 







V 



s z — X 
sw — y 

which evaluated on the periodic orbit T reads for: 
/ -2 

k( 7 (0) = 



\ 

s k — s z (x + z) —l — sz(y + w) 

l — sw(x + z) sk — sw(y + w) / 



\ 

1 — 2scos(t) sin(t) 
s (A;-2sin 2 (t)) / 



(s-l)cos(t) s(k - 2cos 2 (t)) 
\ (s - 1) sin(t) 1 - 2s cos(t) sin(t) 
The fundamental matrix solution of the linear equation v'(i) = k(7(t)) v(t) is: 
• when 2 — 2s + sk ^ 0, 



/ 



v(f) 



(s - 1) e- 2t 
(s-l)e- 2 * 



e- 2t 

e (2-2s+ks)t _ i 

2-2s + ks 

e (2-2s+ks)t _ ]_ 

2-2s + ks 







cos(t) 
sin(t) 



3 (fe-2)st 



cos(i) 





i kst srin(f) 



\ 



3 (fc-2)st 



sin(t) 



cos(t) 



• when 2 — 2s + sk = 0, we put fc = 2(s — l)/s: 



-2t 



v(f) = 







\ 



(s- l)e" 2 *t cos(t) e- 2 'cos(t) -e^-^sin^) 
^ (s-l)e- 2 't sin(t) e- 2 *sin(t) e 2(s - 1) *cos(t) ^ 
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We can compute the eigenvalues of the matrix v(2ir), which by Theorem |2J corre- 
spond to the characteristic multipliers associated to T. These eigenvalues are: 



• when 2 - 2s + sk ^ 0: e" 4 ^, e 2s( - k ~ 2)n and e 2ksn , 

• when k = 2(s - l)/s: e" 4 ^, e~ 47T and e 4 ^" 1 ) . 

Therefore, for instance in the first case (2 — 2s + sk ^ 0) when s k > or s (k — 2) > 
we have that T is unstable and, in the second case (2 — 2s + sk = 0) when s < 1 
then T is Liapunov unstable. 

Example 2. We give an example related to the Mathieu's equation. The 
Mathieu's equation is a particular case of the Hill's equation and it has the form: 



v"(t) + (a + 2qcos(2t))v(t) 



0. 



where a, q are real parameters. See, for instance, the book [HI pp. 121-131] for 
further information about the Hill's equation. 

We consider the following differential system in R 3 : 



x 

y 



-y + zx/2, 
x + zy/2, 
(-2q(x 2 



(9) 



y 2 ) - a) (x 2 + y 2 



1) +z 2 



where a, q G M. This system has the 27r-periodic orbit T := {^(t) \ < t < 2tt} with 
7(t) = (cost, sint, 0). 

We consider the surfaces given by fi(x,y,z) = 0, i = 1,2, with f(x, y,z) = 
(x 2 + y 2 — l,z) . Their intersection gives the periodic orbit T and the crossings 
over it are transversal, as the following computation shows: 



V/i(7(«)) 




2 cos(t) 2 sin(t) 







V/ 2 (7(*)) 







1 


= -2^0 


X( 7 (t)) 




— sin(i) cos(t) 








We have that Df (x) X(x) 
k(x,y,z) :-- 



k(x) f (x), with the following matrix of cofactors: 




x 2 + y 2 



-2q{x 2 



y ) -a 



Therefore, the cofactor matrix over the periodic orbit reads for: 

k « ( » = ( M(t) I 
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where we denote M(t) := — 2gcos(2t) — a. 

By Theorem |2l the stability of T is given by the eigenvalues of v(2tt), where v(i) 
is the fundamental matrix solution of 

v' n (t) v' 12 (t) \ ( 1 W v n {t) v u (t) 
. v' 21 (t) v' 22 {t) ) " " v M W J V Mt) Mt) 

Hence, 

• v'dlt) = M(t)v\i(t) which is the Mathieu's equation with initial conditions 
«ii(0) = 1 and ^(0) = 0, and 

• v i2{t) = M(t) v\2(t) which is the Mathieu's equation with initial conditions 
Ul2 (0) = and v' 12 (0) = 1. 

We have that the system v'(£) = k.(j(t)) v(t) coincides with the characteristic sys- 
tem associated to the Mathieu's equation. Then the stability of the 27r-periodic 
orbit T := |0 < t < 2n} with ^(t) = (cost, sin t, 0) of system Q coincides 

with the stability of the Mathieu's equation, which is studied in and [3 pp. 
128-130]. 

Example 3. Our third example consists on the study of a periodic solution 
related to rigid body dynamics encountered by Steklov [TT] . We consider the motion 
of a rigid body around a fixed point in a uniform gravity field. We denote the weight 
of the body by W and the distance between the center of gravity and the fixed point 
by t. As described in the work jHj, we can consider two frames of reference both 
with the origin at the moving body. The first frame OXYZ is fixed and has axis 
OZ vertical and upward directed. The second frame of reference is moving solidary 
with the body and its axes Ox, Oy and Oz are directed along the major axes of 
inertia for the point O, with corresponding moments of inertia denoted by a, b and 
c. We denote by p, q and r the components of the angular velocity vector of the 
movement and by 71, 72 and 73, the components of the unit vector in the direction 
OZ written in coordinates Oxyz. Steklov considered the case in which the center 
of mass is located on the major axis of inertia, which we assume to be the Ox-axis. 
The Euler-Poisson equations describe the motion of the rigid body: 

(b-c) . (c-a) We . (a-b) Wi 

V = qr, q = — 7 — pr + ^-73, r = pq 72, , 1fV . 

a b b c c (10) 

7i = r 72 - <?73, 72 = Pl3 ~ ?*7i, 73 = qii - Pl2- 
We are going to assume, without loss of generality and following [9,, that: 

b > c, a + b > c, b + c> a, c + a>b, b > a > 2c. (11) 
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The general Euler-Poisson equations exhibit three first integrals, whose expressions 
are described in [2]. We are going to rewrite them for the particular case (jiup. The 
first one corresponds to the projection of the angular momentum onto the vertical: 

Hi := ap7i + bq^ 2 + cr>y 3 . 

The second first integral that we encounter is the geometric property of the vector 
(71,72,73) to be of constant modulus: 

H 2 := 7i + 72 + 73 - 1- 

We also find the full energy (sum of kinetic and potential energies) of the body as 
first integral: 

1/2,2 2\ „™ (a 2 -2ab-2ac + 2bc)W£ 

H 3 := (a-P + bq 2 + cr 2 ) + Wi lx + ± — ^ . 

2 2(o — a) (a — c) 

Straightforward computations show that, if we denote by X the vector defined by 
system (JTHJ), we have ViZ* • X T =0, for % = 1, 2, 3. 

As described in [H], Steklov in [TT] looked for real particular solutions of system 
(fTUj) satisfying the two relations 72 = P2PQ and 73 = /3 3 rp, with /3 2 and /3 3 suitable 
constants, and he found the following particular periodic solution for system (fTUj). 
Let us define the constant [i := ^JWi/a, the dimensionless "time" v := \i (t + t ), 
with to an arbitrary constant, and the following constants and variable: 



a(a — 2c) a(2b — a) ^ 2 b — a 1 / a 

(b — a)(a — c) , (6 — a) (a — c) ' 6 — c' k\j a — c 

Then, the periodic orbit found by Steklov can be written as: 



p(t) := -/^J^— — ^- cn(z;fc), 7l (t) := 1 — cn 2 (z;fc), 

(a — cj a — c 



g(t) := Hi fo^ c \ sn ( z; ^ 2 (*) := sn(^; fc ) cn (^ fc )> 



rf-/) = /'Jt^t dn(z;*), 7s (t) := -J^f^ cn(z; k) dn(z; k), 

(a — c) y (a — c) 

(12) 

where cn(z; fc), sn(z; fc) and dn(z; fc) are elliptic Jacobi functions of variable z and 
module k. 

We recall that the module k must always satisfy < k < 1 and that these elliptic 
Jacobi functions are defined as: cn(z;k) = cos(&m(z; k)) , sn(z; k) = sin(am(;s:; k)) 
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and dn(z; k) = yl — k 2 sin 2 (am(z; k)), where am(z; k) is the Jacobi amplitude 
defined as the inverse function of the elliptic integral of first kind F(u>; k), that is, 
am(z; k) — w if and only if F(u>; k) = z. 

As stated in 0, the minimal positive period T of this Steklov solution with 
respect to the time t is: 

T = AkK(k) 



we ' 

where K(k) is a complete elliptic integral of the first kind. See pfl E] for further 
information about Jacobi elliptic functions and integrals. We only recall that: 

F(z;k) := / Z ; M =, K(k) := F (—; k 
K ' Jo V±~k 2 sin 2 (£) K ' V2' 

and 

cn 2 {z-k) + sn 2 (^;A;) = 1, dn 2 (z; k) + k 2 sn 2 (z;k) = 1. (13) 

Our goal is to study the characteristic multipliers associated to the periodic orbit 
given by Steklov. Since the system has three functionally independent first integrals, 
it is well known that the 6x6 monodromy matrix computed from the variational 
equations associated to the periodic orbit needs to have at least 4 eigenvalues equal 
to 1. Several works target to the aim of obtaining the remaining two characteristic 
multipliers following the classical first variational analysis. In , a numerical study 
of this particular periodic solution is done and a visualization of it is provided. In 
the work j^j, the problem of orbital stability of the Steklov solution is numerically 
examined. In [7], it is proved that if a > 1, then this Steklov periodic orbit is 
orbitally unstable. In this paper we will see how our method allows to reduce the 
computation of these characteristic multipliers to the study of a linear second order 
differential equation. 

We consider the following five hypersurfaces which contain the periodic orbit 
given by Steklov: 

If 2 , 2 2 \ txt „ (a 2 -2ab-2ac + 2bc)W£ 

fi ■= 77 i a P + bq 2 + cr 2 ) + Wt lx + ^ — ^ = 0, 

2 ' 2(b — a) (a — c) 

(a — b)(a — c) 

h := 72 - We(a-2c) Pq = °' 

(a-b)(a-c) 

h '= 73 " Wi(a-2b) Pr = °' 

;= (g - b) p2 + (6-c) r2 _ W£(a - 2b) = 

(a — 2c) a (a — 6) (a — c) 

(a-c) 2 (b-c) 2 Wi(a - 2c) 

h ■= T^TP + Q + 7 7T7 7 = 0. 

(a — 2b) a (a — b)(a — c) 
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The first hypersurface corresponds to one of the described first integrals. The second 
and third hypersurfaces correspond to the ones looked for by Steklov and the last 
two hypersurfaces also contain the periodic orbit and are independent from the 
previous ones. Note that only the first one of the above hypersurfaces is invariant 
by the flow of (fTUJl . When we substitute the expressions of the parameterization 
(fT2|). and using (JT3J), we get that each ft values identically zero for any t. 

The expressions of the other two first integrals in relation with the five polyno- 
mials fi are: 

ap r . , , ac (a(a - 2c) + 26c) p , 

H t = m h + Hh + erf, + ^_ m l_ c ylt U + 



ab{a{a - 2b) + 26c) p 
2{a - 2c){b - c)W£ 



H2 -{we~wk) fl + {j^hj + f2 ) f2 + ((S) +/3 ) /3 + 

/ —ac a5 (a 3 — a 2 b — 2a 2 c + abc + 2ac 2 — 26c 2 ) p 2 

+ 1 — + (14) 

\{b-c)W£ {a-2b) 2 {a-2c){b-c)W£ V ; 

a 2 c 2 ac(4a6 - 3ac - 26c + 2c 2 ) \ 

+ 4(6 - c) 2 W 2 £ 2 h + 4(6 - c) 2 W 2 £ 2 h ) U + 
( -ab a5(a 3 - 2a 2 b + 2ab 2 - a 2 c + abc - 2b 2 c) p 2 

+ \{b-c)W£ + (a-26)(a-2c) 2 (6-c)^ + 

a 2 6 2 ac(2a6 - 3ac - 26c + 2c 2 ) \ 

+ A(b-c) 2 W 2 £ 2 ^ 4{b-c) 2 W 2 £ 2 /5 ' 

where 5 = (a — 6) (a — c)/(W£). We note that over the periodic orbit, each one 
of these first integrals does not conform a unique leaf since its coefficients over the 
hypersurfaces fi = are not constants. This assertion means that the intersection 
of any combination of four of the five hypersurfaces fi = 0, % = 1,2, 3, 4, 5, and the 
set H\ = (equivalently H 2 = 0) contains more points than T. This is the reason 
why we do not directly use these first integrals in the computations. We have com- 
puted the previous expressions because they will lead us to show the existence of 
two characteristic multipliers equal to 1. 

The crossings of the five hypersurfaces fi = 0, i = 1, 2, 3, 4, 5 are normal over 
the periodic orbit because the value of the determinant of the matrix formed by the 
gradients of each of the fi in the first five rows and the vector field in the last row, 
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all of them evaluated over the periodic orbit (|T2j) is: 



det 
+ 
+ 



2(6 - c)W 2 £ 2 (2a(a - 2c) (a - c) 2a(a - c) 



a(a — c) 5 



(6 - a) 2 
q(tf 



+ 



(26-a)(6-a) 



PW 2 + 



2(26 - a)(a -c) (6 -c) ^ 2 , 2(26 + 2c - 3a) (a - c) 2 (b - c) 

a(6 — a) 
8(a-c)(6-c)jy£ 



(26-a)(6-a)W£ 



r{tf + 



a(26 - a) 



7 2 (t)" 



which is positive for any t. This assertion is true because the conditions (jllj) imply 
that all the coefficients are positive except the term (26 + 2c — 3a), which can be 
positive, negative or zero. If it is positive or zero, we already have that det > 0. If 
it is negative, we are going to show that: 



2a(a - 2c) (a - c) 2(3a - 26 - 2c) (a - c) 2 (6 - c) 



(6 -a) 



(2b- a)(b- a)W£ 



r(t) 2 > 



which ensures that det > for any value of t. We note that —1 < dn(^; k) < 1, 
so we consider any value of t for which dn(z; k) 2 is equal to 1 and, using (|T2"j) and 
some computations, we can bound the previous expression by: 

2a(a - 2c) (a - c) 2(3a - 26 - 2c) (a - c) 2 (6 - c) 2 
(6 - a) 2 (2b-a)(b-a)W£ r ^ ~ 

2a(a - 2c)(a - c) 2(3a - 26 - 2c)(a - c) 2 (6 - c) Wi (26 - a)a 



> 



(6 - a) 2 
2a (26 - a) 
(6 -a) 



(2b-a)(b-a)We 



(6 — a) (a — c) s 



> 0. 



The matrix of cofactors k = (&y) associated to the five previous hypersurfaces, 
that is, 



DfX 



( v/1 ^ 




( h \ 






h 


v/ 3 


X = k 


h 






U 


v v/ 5 y 
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reads for: 



k Xj 
hi 

k 2 4 

hi 

ku 

hi 
hi 



for 
— r 



1,2,. ..,5, 



Wt 



>22 



acr 



2{b-c)Wt 

^32 • 



■ o, k 23 

hb ■= ~ 



2ab — ac + 36c) p 



b{a - 2c) 

(a 2 b - 2ab 2 - 2a 2 c + 2a6c + 26 2 c + 2ac 2 - 26c 2 ) r 



2(a-2c)(6-c)iy£ 



ab — 2ac + 36c) p 



■33 



0. 



Wt ° (a-2b)c 
(2a 2 6 - 2a6 2 - a 2 c - 2a6c + 26 2 c + 2ac 2 - 26c 2 ) q 



2(a-26)(6- c)W£ 



> 35 



-a6g 



2(6 -c)W 



&43 

£52 



/C44 

/c 54 



^45 
^55 



o, 

0. 



h 2 
h 3 



-2{b-c)W£r 
ac 

2(b-c)W£q 
ab 



The fundamental matrix of solutions, evaluated in T, of the linear differential system 
of equations 

v'(t) = k( 7 (t))v(f), 



(15) 



with, 



\ 



J 



( \ 

M*) MO ^24 (t) k 25 {t) 

hi(t) h 2 (t) k u (t) k 35 {t) 
A; 42 (t) 
k 53 {t) 

the matrix of cofactors evaluated over the periodic orbit ()12j) . gives the characteristic 
multipliers associated to the periodic orbit encountered by Steklov. We note that 
the first equation of this system is v[(t) = which corresponds to the fact that 
f\ is a first integral of the system and gives us that t>i(t) needs to be a constant. 
Analogously, the other two first integrals give us relations among the values of the 
solutions of system (JT3J). In particular, from (JHj), Hi gives us that the function: 



hi(t) : 
+ 



ap{t) 



Vl {t) + bq{t)v 2 {t) + cr(t)v 3 {t) + 



Wt 

ac (a(a — 2c) + 26c) p(t) 



ab (a(a - 26) + 26c) pit) , , 



2(a-2b)(b-c)W£ w 2(a-2c)(6-c)W^ 
satisfies h[(t) = when considered over the solutions of system (|15|). And the first 
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integral H 2 gives the function: 



/ -ac a5 (a 3 - a 2 b - 2a 2 c + abc + 2ac 2 - 2bc 2 ) p(t) 2 \ 
+ \{b-c)W£ + {a-2b) 2 {a-2c){b-c)Wi ) V *® + 

{ -ah aS(a 3 -2a 2 b + 2ab 2 - a 2 c + abc - 2b 2 c) p(t) 2 \ 
+ \{b-c)W£ + {a-2b){a-2c) 2 {b-c)Wl ) 

where 5 = (a — b)(a — c)/(W£), which is also constant over the solutions of (JT3J). 

We obtain these constant functions over the solutions of (|15|) because each func- 
tion Hi(j(t)), i — 1,2, is for any t. We can then derive Hi(j(t)) with respect to 
t and we deduce the functions hi(t) taking into account that Vi(t) is related with 
V/i( 7 (t)),* = l,2,...,5. 

We first note that the following constants conform a solution of system (|15|): 

ui(t) = 1, v 2 (t) = 0, v 3 (t) = 0, 
2(a - 2b) . . -2(a - 2c) 



a 2 - 2ab - 2ac + 26c w a 2 - 2ab - 2ac + 26c 

The functions hi(t) and h?,{t) take the constant values and 45/ (a 2 — 2ab— 2ac+2bc), 
with 5 = (a — b) (a — c)/(W£), over this solution, respectively Since this solution 
takes a nonzero value on Vi(t), which needs to be constant, we may get the first 
column of the fundamental matrix solution of (|T3j) from the last four columns and 
this particular constant solution. Hence, we only need to take into account the last 
four columns of the fundamental matrix solution of (fTHj) . 

Let us now consider the fundamental matrix v(t) solution of (j!5|) whose initial 
condition is the identity matrix. We define Vij(t) as the function corresponding 
to row i column j of this fundamental matrix of solutions. Using that the three 
functions vi(t), hi(t) and h,2(t) are constants, we can obtain certain relations among 
the rows of this matrix at the value T. 

We define the constants: 

2(2b - a)(b - c)W£ 2^(b-c)W£ 



^fa{b - a) (a - 2c) 3 / 2 ' ' ' (b - a)Ja~^Tc{a 2 - 2ab + 26c) ' 
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and we get that v(T) equals 



/ 1 \ 

V 2 a(T) v 25 (T) 

ViijT) - 1 f 55 (T) - 1 

Pi P2 

'v 55 (T) - 1 



U2l(T) 


V22(T) 


V2 3 (T) 




V 3 2 (T) 


v 33 (T) 


p\v 3 \{T) 


P\V 32 {T) 


Pl (v 33 (T)-l) 


P2V 3 l(T) 


P2V 3 2(T) 


p2 (v 33 (T) - 1) 



V M (T) pi 

vu(T)-i 



P2 



The computation of the characteristic polynomial associated to this matrix gives: 
det (v(T) -XI) = -(A - l) 3 (A 2 + B A + C) , 



where 
and 



B := 2 - V22{T) - v 33 {T) - v iA {T) - v 55 (T), 

C := v 2 2{T) K(T) + v u (T) + v 55 (T) - 2) + 
-V32 {T) (v 23 (T) + Piv 2 a{T) + P2V2b{T)). 

We note that each one of the first integrals gives an eigenvalue equal to one. More- 
over, using Liouville's formula, we know that the product of all the eigenvalues of 
v(T) is equal to one because the trace of matrix k(7(t)) is identically zero and v(t) 
is a fundamental solution of (|15p. Hence, we get that C — 1. Then, we have that 
the roots ofA 2 + £?A + l = are: 

• if B 2 > 4, then the roots are both real, and one of modulus greater than 1 
and the other with modulus lower than 1: so the Steklov solution is unstable, 

• if B 2 < 4, then the roots are both of modulus equal to one and the charac- 
teristic multipliers do not decide if the Steklov solution is unstable or not. 

Let us go back to the differential linear system We know that vi(t) is 

constant, and we denote its value by v w := i>i(i). Analogously, we denote by 
^10 := ^i(^) and h 2 o := h 2 {t), the values of the other two constants. We can 
determine v 2 {t) and v 3 (t) from the last two equations of (fl3j) and also from the two 
relations hi(t) = h w and ti2(t) = /i 2 o- We equate the two expressions of V2(t) and 
v 3 (t) obtained in these two different ways and we get a linear differential system of 
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equations for v±{t) and v 5 (t) which reads for: 

ddn(z;k) W£(2b - a)dn(z; k) 

■v 4 (t) - 



dt 
+ [b- 

+ 



a(a — c)(b — c) 
(26 — a) (ab — ac — be) 



ab — 2ac + 26c . . 2 , 
c 2b -a cn ^ k ) )M t ) + 



(a — 2c)c 
2(6 - c)y/We(b-a)(a-2c) 



ac\j2b — a 



cn(z; kyj v 5 (t) + 

/TXT 2(6 — c) 2 
cn(z; ft) ftio H Uio cn(z; ft) + 



+ — (tW/l 20 - 2uio) 
a 



0, 



<9sn(z;ft) , ^ Wlia — 2c)sn(z; ft) 



9* 
+ (6 + 

+ 



a(b — a) (6 — c) 
(6 — a)(2ab — ac — 26c) 



c - ro(n (26- c r c) cnfet)2| " 4(f)+ 



(a — 2c) (a — c) 
-2(6 - c) v /W£(b-a)(2b-a 



aby/a — 2c 

+ — (W£h 20 - 2v w ) 
a 



cn{z; kfj v 5 (t) + 
cn(z; ft) h w 



2 (6-a)(6-c) _,_.^ 2 



6(a — c) 



cn(z; ft) uio + 



0, 



where r$ = (6 — a)(a — 2c)/(6(a — c)). Equating fs(t) from the first equation and 
substituting its value in the second one, we get a second order linear differential 
equation for v±{t) whose fundamental set of solutions would let us to the compu- 
tation of the characteristic multipliers associated to the Steklov periodic orbit. We 
define 

(26 - a)S 1 



u(t) := 6 



— cn(^; A;) 2 , Si := ab — ac — bc 



(a — 2c)c 

and this second order linear differential equation for v^(t) is: 
V(t) 2 ddn(z-k) 



v'l(t) - 

where 
A (t) 

we 

A nh (t) :-■ 



v' A (t) + A (t)v 4 (t) + A nh (t) = 0, (16) 



u(t) dn(z; k) dt 

2(b-c)(a 2 -Ab(a-c)) dn(z;k) 2 4 a 2 - 2(5 1 + 2ac) 



a(a — 2c) (a — c) 
-2(26 -a) 2 



+ - + 

a 



cn(z; ft) 2 , 



a5(a — 2c)c 2 uj(t) 
with 5 := (a - 6) (a - c)/(W£) and 



u(t) a bc(a — c) 

«3 (t) hio + a 4 (t) v io + a 5 (t) (2v 10 - W£ h 2 o) ) , 



a 3 (t) : = 



W£(b - a) 



(26 -a) (a -2c) 



[a - 2c)(3ac(6 + c - a) + 56c(6 - c)) 
(26 - a) 



+ 
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+2ab{2c-a) + 



5 1 ((b-a) 2 + 5 1 ) 



cn(z; k) 2 I cn(z; k) 



a 4 (t) := -a(a - 2c) f — ^— + (6 - 2c) cn(z; A;) 2 + % — cn{z; fc)- 

\2b — a b[a — 2c) 

, . abc c 2 (a 2 — ab — b 2 — ac+ be) (a — 2b + 2c) Si , , . 9 

a 5 (t) := — + ^ '--I '—cn.{z;kf. 

w 2 2b -a 2 K 1 J 



In short, we have reduced the problem of computing the characteristic multipli- 
ers for the periodic orbit (JT2jl to the study of the second-order linear differential 
equation (|16)1. The classical approach to this problem is to compute the characteris- 
tic multipliers via the first variational equations, thus involving a linear differential 
system of order 6. Our method starts with a linear differential system of order 5, see 
(I15|) . In this particular problem, we know three first integrals which let us reduce 
the order by 3, getting the second-order linear differential equation (fTB|l. We note 
that the use of the first integrals is not trivial, since we need to write them in terms 
of the considered hypersurfaces, see (1 14 J) , and then relate them to the variables of 
the linear differential system (fTK|) . 
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